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ABSTRACT 

We  present  the  results  of  differential  phase  experiments  done  with  data  from  the  Navy  Prototype  Optical  In¬ 
terferometer  (NPOI).  We  take  advantage  of  the  fact  that  this  instrument  simultaneously  records  16  spectral 
channels  in  the  wavelength  range  550-850nm,  for  multiple  baselines.  We  discuss  the  corrections  applied  to  the 
data,  and  show  the  results  obtained  for  Vega  and  the  Be  star  (3  Lyrae. 
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1.  INTRODUCTION 

Differential  phase  is  a  powerful  technique  used  to  measure  small  phase  variations  as  a  function  of  wavelength,  due 
to  sources  with  chromatic,  asymmetric  structures.  Previous  applications  of  differential  phase  include  the  study 
of  Ha  disks  and  envelopes  around  stars,12, 13  and  the  detection  of  phase  fluctuations  due  to  a  binary  companion.1 
Since  this  technique  relies  on  simultaneous  multiwavelength  phase  measurements,  it  is  ideally  suited  for  NPOI 
observations.  In  the  following  we  describe  the  first  tests  of  this  technique  on  NPOI  data  from  a  star  with  a 
circumnuclear  Ha  disk  and  from  an  asymmetric  star. 

2.  OBSERVATIONS  AND  DETERMINATION  OF  DIFFERENTIAL  PHASES 

The  observations  presented  in  this  paper  were  obtained  with  the  NPOI2  on  the  nights  of  2004-09-01  and  2005-05- 
19.  These  observations  were  done  with  2  spectrographs,  simultaneously  recording  16  channels  in  the  wavelength 
range  550-850nm,  and  3  baselines  per  spectrograph.  The  first  set  of  observations  used  baselines  with  lengths 
between  18  m  and  64  m,  while  in  the  second  run  the  longest  baseline  was  of  the  order  of  ~53  m.  The  night  of 
2004-09-01  was  used  to  observe  several  scans  of  Vega  (a  Lyrae),  its  calibrators  (7  Lyrae)  and  other  stars,  while 
on  2005-05-19  we  obtained  9  scans  of  (3  Lyrae  and  its  calibrator  7  Lyrae.  Each  scan  had  a  duration  of  30  s,  and 
the  observations  of  the  program  star  were  interleaved  with  observation  of  its  calibrator. 

The  data  reduction  started  with  the  coherent  integration  of  complex  visibilities,  which  was  done  applying 
the  algorithm  developed  by  Ref.  5  (see  Refs.  6,7  for  a  new  coherent  integration  algorithm).  This  technique  uses 
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Figure  1.  This  figure  presents  baseline  phases  (AEO-ACO,  18  m)  of  7  Lyrae  as  a  function  of  wavenumber.  The  left 
panel  shows  the  raw  phases  from  the  first  scan  observed  in  the  night  of  2005-05-19.  Each  individual  line  corresponds  to 
a  coherently  integrated  subscan  with  a  duration  of  200  ms.  The  right  panel  shows,  as  a  white  line,  the  average  of  all 
calibrator  scans  obtained  on  this  night,  including  7  Lyrae.  The  black  lines  correspond  to  all  200  ms  coherently  integrated 
subscans  of  7  Lyrae  observed  on  that  night  (1728  in  total),  corrected  for  differential  atmospheric  delays  to  match  the 
white  line. 

the  deviation  of  the  calculated  fringe  position  from  the  estimated  geometrical  value  to  rotate  the  phasors  onto  a 
common  fringe.  The  aligned  phasors  are  then  averaged  into  200  ms  subscans  (100x2  ms  instrumental  integration 
times),  and  the  phases  are  calculated  for  the  individual  channels. 

The  results  of  this  procedure  are  presented  in  the  left  panel  of  Fig.  1,  where  we  show  the  phases  for  each 
200  ms  subscan  of  the  first  scan  of  the  calibrator  star  7  Lyrae.  We  can  see  that  the  individual  subscan  phases 
follow  smooth  curves  with  similar  shapes.  Another  important  result  that  can  be  seen  in  this  figure  is  the  fact 
that  the  subscans  with  high  phase  values  in  the  red  have  lower  ones  in  the  blue,  and  vice-versa.  Since  this  star  is 
not  resolved  by  these  observations  it  should  have  zero  phase  on  average,  so  the  observed  phases  can  be  attributed 
to  uncompensated  glass  and  differential  airpath  between  the  two  beams.  The  effect  of  uncompensated  glass  is 
not  expected  to  vary  significantly  over  the  course  of  the  night,  or  even  over  the  course  of  several  days,  so  we 
can  identify  this  effect  as  being  responsible  for  the  overall  shape  of  the  curves  in  the  left  panel.  The  variable 
component,  responsible  for  the  subscan  to  subscan  variations  can  be  identified  with  the  effect  due  to  differential 
airpath  delay  between  the  two  stations.  This  effect  adds  a  quadratic  phase  as  a  function  of  wavenumber  to  the 
data. 

In  the  right  panel  of  Fig.  1  we  show  that  fitting  and  subtracting  a  curve  of  the  form  y  =  a  +  bx  +  cx 2  ( x  is 
wavenumber),  to  the  variable  phase  component  of  the  subscans,  we  can  reduce  them  to  a  common  curve,  shown  in 
white.  This  curve  corresponds  to  the  average  of  all  calibrator  scans  observed  throughout  the  night.  By  averaging 
all  these  scans  we  basically  eliminate  the  differential  airpath  delay  effects,  under  the  reasonable  assumption  that 
the  mean  airpath  difference  over  the  course  of  a  night  is  zero,  and  are  left  only  with  the  instrumental  phases 
(uncompensated  glass).  The  black  lines  show  all  the  200  ms  subscans  of  7  Lyrae  observed  on  that  night  (1728  in 
total),  after  fitting  and  subtracting  the  function  described  above.  Notice  that  all  lines  cluster  around  the  white 
line. 
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Figure  2.  This  figure  shows  the  phases  of  the  Be  star  j3  Lyrae  on  the  long  baseline  E06-AC0  (34.2  m),  observed  on 
2005-05-19.  The  top  left  panel  shows  the  instrumental  phases  (dashed  line)  and  the  average  phases  of  9  individual  scans 
(solid  lines),  as  a  function  of  wavenumber.  The  vertical  doted  line  shows  the  location  of  Ha.  The  top  right  panel  shows 
the  phases  of  each  scan  after  the  subtraction  of  the  instrumental  contribution.  The  bottom  panel  shows  the  residuals 
obtained  after  fitting  a  curve  of  the  form  y  =  a  +  bx  +  cx 2  to  the  continuum  points  of  each  individual  curve  in  the  top 
right  panel,  and  subtracting  it  from  all  the  channels,  including  Ha.  The  uncertainty  in  the  phases  is  of  the  order  of  0.5°. 
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2.1.  (3  Lyrae 

The  results  presented  in  Fig.  1  show  that  in  the  case  of  calibrator  stars,  which  are  symmetric  and  the  phase  is 
intrinsically  zero,  we  can  determine  and  correct  the  instrumental  phases,  as  well  as  the  effects  due  to  differential 
airpath  between  two  stations.  We  show  in  Fig.  2  the  application  of  this  technique  to  one  of  the  baselines  of 
the  Be  star  f3  Lyrae,  where  we  aim  at  detecting  the  phases  of  the  Ha  channel.  This  is  a  binary  system,4  only 
partially  resolved  by  our  observations,  where  both  components  have  similar  magnitudes  and  colors,  so  the  system 
is  approximately  symmetric  and  should  have  phases  close  to  zero  in  the  continuum  channels.  However,  in  the 
case  of  the  Ha  channel  we  expect  to  see  a  phase  signature.  This  line  originates  in  a  disk  ~2  mas  in  diameter, 
that  surrounds  only  one  of  the  stars  in  the  system.  As  a  result  of  this  geometry,  the  Ha  photocenter  is  displaced 
relative  to  the  photocenter  of  the  binary  system,  and  should  show  a  differential  phase  relative  to  the  continuum. 

We  followed  the  same  reduction  steps  discussed  above,  but  analyze  the  average  of  each  30  s  scan,  instead  of 
the  individual  200  ms  subscans.  We  take  advantage  of  the  fact  that  the  differential  airpath  effects  are  additive, 
so  we  can  fit  and  subtract  the  effective  contribution  of  this  effect  to  the  average  scan  of  a  source.  The  top  left 
panel  of  Fig.  2  shows  the  average  phases  of  the  9  scans  of  /3  Lyrae,  the  instrumental  phase  obtained  from  the 
average  of  all  calibrator  scans  and  also  indicate  the  location  of  the  Ha  line.  The  top  right  panel  shows  the  phases 
obtained  after  the  subtraction  of  the  instrumental  phases,  which  leaves  only  the  phases  due  to  the  source  and 
differential  airpath.  We  can  see  that  the  phases,  as  in  the  case  of  the  calibrator  star  7  Lyrae,  vary  smoothly  with 
wavelength,  except  for  a  clear  deviation  in  the  Ha  channel.  The  correction  for  the  effect  of  differential  airpath  is 
done  by  assuming  that  the  phases  should  be  zero  in  the  continuum  channels,  so  their  phases  were  fitted  as  in  the 
case  of  7  Lyrae,  and  subtracted  from  all  channels,  including  Ha.  The  residuals  from  this  operation  are  presented 
in  the  bottom  right  panel  of  the  figure,  where  we  can  see  that  the  final  continuum  phases  cluster  around  zero 
and  the  Ha  phases  show  a  clear  signal. 

The  differential  phases  of  (3  Lyrae,  for  all  baselines  observed  in  2005-05-19,  are  presented  in  Fig.  3.  As 
in  Fig.  2,  we  can  see  that  the  strongest  phase  signal  is  in  Ha,  with  the  continuum  differential  phases  usually 
clustering  around  zero.  We  also  see  some  deviations  from  zero  around  l.l/iiiW  and  1.7/zm-1,  which  correspond 
to  the  neutral  Helium  lines  at  706.5  nm  and  587.6  nm,  respectively.  Other  deviations  from  zero  in  the  continuum 
channels  can  be  due  to  atmosphere  or  instrument  effects,  and  require  further  investigation.  A  particularly 
interesting  result  from  this  figure  is  the  fact  that  the  larger  phase  variations  occur  along  the  E-W  baselines.  This 
indicates  that  the  binary  and  the  disk  are  extended  along  this  direction,  consistent  with  the  non-resolution  of 
the  binary  by  the  GI2T  interferometer,  which  had  a  N-S  orientation. 

2.2.  VEGA 

Another  application  of  the  differential  phases  technique  is  the  study  of  stellar  structure.  Here  we  present  the 
particularly  interesting  case  of  Vega  (a  Lyrae),  the  primary  standard  calibration  star.  Recent  observations  by 
NPOI10  and  CHARA3  showed  that  this  is  a  rapidly  rotating  star  (~93%  of  the  breakup  velocity),  seen  almost 
pole-on  ( i  =  4.54°  along  p.a.  =  8.6°).  The  effects  of  the  rotation  cause  the  star  to  have  an  equatorial  radius 
1.246  times  larger  than  the  polar  one,  and  a  temperature  drop  of  ~2400  K  from  the  poles  to  the  equator.  These 
properties  make  the  star  look  asymmetric  when  seen  projected  in  the  sky,  an  effect  that  produces  a  strong  phase 
signature  at  wavelengths  around  the  first  null  crossing  (this  signature  was  used  by  Ref.  10  to  fit  the  properties 
of  this  star) . 

We  show  in  the  top  panels  of  Fig.  4  and  Fig.  5  the  predicted  differential  phases  of  Vega  for  the  baseline 
W07-AE0  (64.4  m),  observed  on  2004-09-01  at  hour  angles  of  2.92  h  and  3.55  h,  respectively.  This  figure  shows 
the  first  null  crossing  around  700  nm  and  670  nm  (phase  flip  of  180°),  where  we  can  also  see  the  large  phase 
deviations  due  to  the  asymmetry  of  the  star.  The  bottom  panels  of  the  figures  show  the  observed  differential 
phases,  after  a  first  order  subtraction  of  the  differential  airpath.  We  can  see  that  we  detect,  in  both  cases,  a 
phase  signature  similar  to  the  predicted  one,  but  with  some  small  differences.  In  the  case  of  Fig.  4,  the  observed 
phase  flip  happens  at  a  higher  wavelength  relative  to  the  predicted  value.  Also,  the  amplitude  of  the  phase 
flips  does  not  match  the  predicted  values  in  both  cases.  We  attribute  these  differences  to  wavelength  calibration 
errors  and  cross  talk  effects,11  and  are  currently  investigating  other  possible  explanations. 
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Figure  3.  Residual  phases  of  f5  Lyrae,  calculated  for  6  different  baselines  observed  on  2005-05-19  (E06-AC0=34.2  m, 
AW0-E06=53.2  m,  AW0-AC0=22.2  m,  AE0-AC0=18.9  m,  AN0-AC0=22.9  m,  AE0-AN0=34.9  m)  .  Each  row  corresponds 
to  a  different  spectrograph.  Each  solid  line  is  a  different  scan  and  the  dotted  lines  indicate  the  position  of  Ha. 


3.  SUMMARY 

We  presented  the  results  of  differential  phase  experiments  done  with  multiwavelength  data  from  the  NPOI.  We 
discussed  the  reductions  applied  to  the  data,  the  correction  of  phases  due  to  instrumental  effects  and  differential 
airpath  between  stations.  We  presented  the  application  of  this  technique  to  the  Be  star  (3  Lyrae,  which  shows  a 
strong  differential  phase  signature  in  the  Her  channel,  particularly  along  the  E-W  baselines,  consisted  with  pre¬ 
vious  results.  Further  applications  of  this  technique  to  Be  stars  include  the  determination  of  complex  visibilities 
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Figure  4.  The  top  panel  shows  the  predicted  differential  phases  of  Vega,  based  on  NPOI  results  from  Ref.  10  (inclination 
and  position  angle  of  the  pole  are  indicated  on  the  top  of  the  panel),  for  observations  with  the  baseline  W7-AE  (64.4  m) 
on  the  night  of  2004-09-01  at  an  hour  angle  of  2.92  li.  The  bottom  panel  shows  the  corresponding  observed  differential 
phases,  after  the  correction  for  instrumental  effects  and  a  first  order  subtraction  of  the  differential  airpath  effects. 


and  the  imaging  of  their  Her  disks  using  standard  radio  interferometry  software,  presented  in  this  conference  by 
Ref.  9.  We  also  presented  the  case  of  Vega,  a  rapidly  rotating  star  seen  almost  pole-on.  We  show  that  we  can  see 
the  phase  signature  due  to  the  asymmetric  structure  of  the  star,  but  we  still  need  to  understand  some  instrument 
issues,  like  cross  talk,  in  order  to  further  improve  the  use  of  this  technique  to  other  stars.  This  technique  will  be 
useful  for  the  detection  of  asymmetric  structures,  due  to  rotation,  flares  or  spots,  in  other  stars. 
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Figure  5.  Same  as  Figure  4  for  hour  angle  3.55  h. 
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